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The study of good nonregular fractional factorial designs has re- 
ceived significant attention over the last two decades. Recent research 
indicates that designs constructed from quaternary codes (QC) are 
I very promising in this regard. The present paper shows how a trigono- 

metric approach can facilitate a systematic understanding of such QC 
designs and lead to new theoretical results covering hitherto unex- 
plored situations. We focus attention on one-eighth and one-sixteenth 
fractions of two-level factorials and show that optimal QC designs 
often have larger generalized resolution and projectivity than com- 
^ , parable regular designs. Moreover, some of these designs are found 

00 ' to have maximum projectivity among all designs. 

a^ 

CN ' 1. Introduction and preliminaries. Fractional factorial designs play a key 

lO '. role in efficient and economic experimentation with multiple factors and have 

' gained immense popularity in various fields of application such as enginee- 

ring, agriculture and medicine. These designs are broadly categorized as 
regular and nonregular depending on whether or not they can be generated 
via defining relations among the factors. In regular designs, any two factorial 
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effects are either mutually orthogonal or completely aliased, and the criterion 
of maximum resolution [2] and its refinement, minimum aberration (MA) 
[10], are commonly used in discriminating amongst these designs. We refer 
to [14, 22] for detailed surveys and extensive references on regular designs. 

The last two decades, especially the last ten years, have witnessed a sig- 
nificant spurt in research on nonregular designs. The case of two-level factors 
has received particular attention. The notions of resolution and aberration 
have been generalized, with statistical justifications, to these designs; see 
[7, 9, 13, 19, 20, 27, 28]. More recent work on nonregular designs and re- 
lated topics include [4] giving theoretical results on generalized MA designs, 
[12] giving a catalog of generalized MA designs, [24] on moment aberration 
projection designs, [23] on designs obtained from the Nordstrom and Robin- 
son code, and [17] on a complete classification of certain two-level orthogo- 
nal arrays. It is well recognized now that although nonregular designs have 
a complex aliasing structure, they can outperform their regular counterparts 
with regard to resolution or projectivity, and this is one of the principal mo- 
tivating forces for the current surge of interest in these designs. For more 
details, see [25] giving a state-of-the-art review of nonregular designs with 
a comprehensive list of references. 

A recent major development in nonregular two-level designs has been the 
use of quaternary codes (QC) for their efficient construction. The resulting 
two-level designs are hereafter called QC designs. While QCs are known to 
yield good binary codes in coding theory [11], QC designs have been seen to 
be attractive with regard to resolution, aberration and projectivity. More- 
over, as noted in [25], these designs are relatively straightforward to con- 
struct and have simple design representation. Xu and Wong [26] pioneered 
research on QC designs and reported theoretical as well as computational 
results. Phoa and Xu [16] obtained comprehensive analytical results on quar- 
ter fraction QC designs and showed that they often have larger resolution 
and projectivity than regular designs of the same size. 

The present paper aims at extending [16] to more highly fractionated 
settings. A serious hurdle in this regard is that the inductive proofs in [16] 
become unmanageable when one attempts to go beyond quarter fractions. 
A trigonometric representation for QC designs is employed here in order to 
overcome this difficulty. This approach is found to be quite convenient for 
one-eighth and one-sixteenth fractions which form our main focus. Earlier, 
Phoa [15] in his unpublished Ph.D. dissertation reported partial results on 
one-sixteenth fraction QC designs under a certain assumption. An advantage 
of our approach is that it involves no such restrictive assumption and enables 
us to obtain unified and more comprehensive results on these fractions. As 
discussed in the concluding remarks, the trigonometric formulation holds 
the promise of being applicable to even more general settings as well. 
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In Section 2, we introduce the trigonometric approach and present results 
pertaining to generahzed resolution and wordlength pattern (WLP) for the 
case where the number of runs, N, is an even power of 2. The corresponding 
results when N equals an odd power of 2 appear in Section 3. Section 4 
dwells on projectivity and some directions for future work are indicated 
in Section 5. Satisfyingly, with one-eighth and one-sixteenth fractions, at 
least over the range covered by our tables, (a) the same design turns out 
to be optimal among all QC designs with respect to all the commonly used 
criteria like resolution, aberration and projectivity, and (b) such an optimal 
design is often seen to have higher resolution and projectivity than what 
regular designs can achieve. The point in (b) reinforces the findings in [16] 
for quarter fractions but that in (a) is in contrast to what they observed 
in their setup. It is also seen that some of our optimal QC designs have 
maximum projectivity among all designs. 

Before concluding this section, we reproduce some definitions from [16] 
for ease in reference. A two-level design D with N runs and q factors is 
represented by an x g matrix with entries ±1, where the rows and columns 
are identified with the runs and factors, respectively. For any subset S = 
{ci, . . . , Cfc} of k columns of D, define 



where Csj is the sth entry of column Cj. The Jk{S;D) values are called the 
J- characteristics of design D; cf. [9, 20]. Following [8], the aliasing index 
of S is defined as Pk{S;D) = \ Jk{S; D)\/N . Clearly, < Pk{S;D) < 1. If 
Pk{S;D) = 1, then the columns in S are fully aliased with one another and 
form a complete word of length k and aliasing index 1. If < pk{S;D) < 1, 
then these columns are partially aliased with one another and form a partial 
word of length k and aliasing index pk{S]D). Finally, if pk{S;D) = 0, then 
these columns do not form a word. 

Let r be the smallest integer such that inax^s=r Pr{S; D) > 0, where # 
denotes cardinality of a set and the maximum is over all subsets S" of r 
columns of D. The generalized resolution [9] of D is defined as 



N 





s=l 





For 1 < < g, let 



(1.3) 




#S=k 



The vector {Ai(D), . . . , Aq(D)) is called the generalized WLP of D. The 
generalized MA criterion [27], also known as minimum G2 aberration [20], 
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calls for sequential minimization of {Ai{D), . . . ,Aq{D)). When restricted to 
regular designs, generalized resolution, generalized WLP and generalized 
MA reduce to the traditional resolution, WLP and MA, respectively. For 
simplicity, we use the terminology resolution, WLP and MA for both regular 
and nonregular designs. 

Following [1], the design D is said to have projectivity p if every p-factor 
projection contains a complete 2^ factorial design, possibly with some points 
replicated. Evidently, a regular design of resolution R has projectivity R — 1. 
As shown in [9], a design with resolution R> r has projectivity greater 
than r. 

2. Quaternary code designs in 2^*^ runs. 

2.1. One-sixteenth fractions. In the spirit of [16], let C be the QC given 
by the n x (n + 2) generator matrix [u v In], where u = {ui, . . . ,Un)' and 
V = {vi, . . . jVnY are n x 1 vectors over Z4 = {0,1,2,3} (mod 4), /„ is the 
identity matrix of order n over Z4, and the primes stand for transpose. 
The code C, consisting of 4"(=2^"') codewords, each of size n + 2, can be 
described as 

(2.1) C = {{a'u, a'v, ai, . . . , a„) : oi, . . . , € Z4}, 

where a = (oi, . . . , a„)', and a'u and a'v are reduced mod 4. The Gray map, 
which replaces each element of Z4 with a pair of two symbols, transforms C 
into a binary code D, called the binary image of C. For convenience, the two 
symbols are taken as 1 and —1, instead of the more conventional and 1. 
Then the Gray map is defined as 

(2.2) 0^(1,1), 1^(1,-1), 2^(-l,-l), 3^(-l,l). 

With its codewords as rows, D is a 2^'" x (2?! + 4) matrix having entries ±1. 
Indeed, with columns and rows identified with factors and runs, respectively, 
D represents a design involving 2n + 4 two-level factors and 2^" runs. In this 
sense, D will be referred to as a 2(^"'+^)~^ QC design. 

A representation of D using trigonometric functions facilitates the study 
of its statistical properties which depend on the choice of u and v. Since the 
pair (^/2sin(f + ffe),\/2cos(f + f A;)) equals (1,1), (1,-1), (-1,-1) and 
(-1,1) for A; = 0,1, 2 and 3 (mod4), respectively, by (2.1) and (2.2), the 2^"" 
runs in D can be expressed as 
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(2.3) 



TT vr 
sinl 4 + 2"^ 



, cos 



TT TT 
4 + 2"" 



ai, . . . ,a„ G Z4. 



Denote the 2n + 4 factors in D by Fi, . . . , i<4, Fn, F12, . . . , -Fni; -Pn25 in con- 
formity with the ordering in (2.3), that is, V^sm{j + f a'u) and \/2cos(^ + 
^a'u) are the levels of Fi and F2, and so on. 

The factors Fi,...,F4, with levels dictated by u or v, require special 
attention. Prom this perspective, for any nonempty collection of factors (or 
equivalently, columns of D), let = 1 or according as whether F/^. is 
included in the collection or not, 1 < /c < 4. With this notation, the collection 
is said to be of the type x = xiX2X^X4,. Thus, for any binary 4-tuple x = 
xiX2X^Xi, a typical collection of type x consists of factors F^ with x^ = 1 
(1 < A; < 4), and also factors Fji,Fj2 {j £ Si),Fj2 (j G 5*2) and Fji (j G ^3), 
where Si, 82,83 are any disjoint subsets of {1, . . . , n}. The total number of 
factors in the collection is then m + X , where X = xi + X2 + X3 + X4,m = 
2ni + 77.2 + ^3 and nj = ^8j (1 < j < 3). Here 81,82,83 can be empty sets 
as well but if 2; = 0000 then at least one of them is nonempty, for otherwise, 
the collection contains no factor at all. Prom (1.1), (2.3) and the definition 
of aliasing index, it follows that the aliasing index of a collection of type x 
{=xiX2X3Xi) as described above is given by |y(x)|, where 



(2.4) 
with 



3 3 
ai=0 a„=0 



(2.5) 



and 



(2.6) 



[x;ai,...,an) 



X sin-^^ ( - + —CLV I cos"""* ( ^ + -^a'"" 



X ii{ai, ...,an) 



ip{ai,...,an) 



TT < sinf — I — Gj I cosf — I — a 



TT vr 
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By (2.4)-(2.6), for any fixed x, the value of V{x) depends on the sets 
51,52,53 in addition to u and v. Any choice of 5i,52,53 that makes V{x) 
nonzero entails a word of length 171 + X and aliasing index Such a 

word will be called a word of type x. 

We now present Theorem 1 below giving an account of words of all possible 
types. For x = 0101, this result has been proved in the Appendix. The proofs 
for all other x are similar and occasionally simpler. In particular, the case 
X = 0000 is evident from the presence of In in the generator matrix [u v In] 
of C. Some more notations will help. With reference to the vectors u and v, 
let 





A 


= {3--^<J<n,Uj =k,Vj = s}, 


(2.7) 






fks 


= #Aks, 0<k,s<3, 




Ai 


= /lO + /305 A2 = /oi + /o3; A3 = /i2 + /32, 


(2.8) 


A4 


= /21+/235 A5 = /ii + /33, A6 = /l3 + /31> 




At 


= /02, As = /205 A9 = /22, Aio = /oO; 




h 


= 2(A4 + As + A9) + Ai + A3 + As + Ae, 




h 


= 2(A3 + At + A9) + A2 + A4 + As + Ae, 




h 


= 2(A2 + As + A9) + Ai + A3 + As + Ae, 




k 


= 2(Ai + At + A9) + A2 + A4 + As + Ae, 


(2.9) 






h 


= 2(Ai + A3 + As + Ae), k = 2(A2 + A4 + As + Ae), 




h 


= 2(Ai + A2 + A3 + A4), h = 2(At + As) + Ai + A2 + A3 + A4 




h 


= 2(A5 + At + As) + Ai + A2 + A3 + A4, 




ho 


= 2(A6 + At + As) + Ai + A2 + A3 + A4, 




Pi 


^ -^y2<(V2)(Ai+A3+A5+A6))^ ^ ]^/2<(^/^)(^2+A4+A5+A6)>^ 


(2.10) 


6 


^ ^/2<(V2)(Ai+A3)>^ ^2 ^ ^/2<(V2){A2+A4))^ 






^ ;^^2<(V2)(Al+A2+A3+A4 + l)>^ 



where (y) is the largest integer not exceeding y. By (2.7), fks equals the 
frequency with which (ks) occurs as a row of the n x 2 matrix [u v]. Also, 
each of the quantities introduced in (2.8)-(2.10) is uniquely determined by 
these frequencies. 



Theorem 1. With reference to the 2(2"+^) QC design D, the following 
hold: 

(a) For X = 0000, there is no word of type x. 



QUATERNARY CODE DESIGNS 



7 



(b) For each o/ x = 0100 and 1000, there are l/pf words of type x; every 
such word has aliasing index pi and length /i + 1. 

(c) For each of x = 0001 and 0010, there are l/p2 words of type x; every 
such word has aliasing index p2 and length I2 + I. 

(d) For each 0/ 2; = 0111 and 1011, there are l/pf words of type x; every 
such word has aliasing index pi and length /s + 3. 

(e) For each of x = 1101 and 1110, there are I/P2 words of type x; every 
such word has aliasing index p2 and length Z4 + 3. 

(f) For X = 1100, there is one word of type x, with aliasing index 1 and 
length 1^ + 2. 

(g) For X = 0011, there is one word of type x, with aliasing index 1 and 
length /g + 2 . 

(h) For X = 1111, there is one word of type x, with aliasing index 1 and 
length Ij + 4. 

(i) For each of x = 0101 and 1010, 

(11) if X5 + Ae = 0, then there are 1/(^1^2) words of type x; every such 
word has aliasing index ^1^2 o-n-d length /§ + 2; 

(12) i/ A5 + Ag > and Ai + A2 + A3 + A4 = 0, then there is one word of 
type x, with aliasing index 1 and length Iiq + 2; 

(13) i/ A5 + Ae > and Ai + A2 + A3 + A4 > 0, then there are 1/^^ words 
of type X, each with aliasing index S,; half of these words have length 
Ig + 2 and the rest have length /lo + 2. 

(j) For each of x = 0110 and 1001, 

(jl) i/A5 + A6 = 0, then the same conclusion as in (il) holds; 

(j2) i/ A5 + Ag > and Ai + A2 + A3 + A4 = 0, then there is one word of 

type x, with aliasing index 1 and length I9 + 2; 
(j3) i/ A5 + As > and Ai + A2 + A3 + A4 > 0, then the same conclusion 

as in (i3) holds. 

Since ^ = 1 when Ai + A2 + A3 + A4 = 0, one can merge (i2), (i3), (j2) and 
(j3) above when words of all types are considered together. Summarizing 
Theorem 1, we thus get the next result. 

Theorem 2. With reference to the 2(2"+4)-4 q(j design D, the following 
hold: 

(a) There are 4//?f words each with aliasing index pi, half of these words 
have length li + l and the rest have length ^3 + 3. 

(b) There are 4/ p^ words each with aliasing index P2; half of these words 
have length ^2 + 1 oind the rest have length ^4 + 3. 

(c) There are three words each with aliasing index 1; these have lengths 
k + 2,1^ + 2 and I7 + A. 

(d) In addition, 
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(dl) if A5 + Ae = 0, then there are 4/ (^^^2) words each with aliasing 

index ^1^2 and length Ig + 2; 
(d2) if X5 + As > 0, then there are 4/^^ words each with aliasing index 

^; half of these words have length I9 + 2 and the rest have length 

ho + 2. 

Theorem 2, in conjunction with (2.9) and (2.10), shows that the resolution 
and WLP of the design D depend on u and v only through Ai,...,Aio. 
Indeed, for any given u and v, Theorem 2 readily yields these features of D. 
This is illustrated below. 

Example 1. With n = 3, let n = (2,1,1)' and v = (1,1,3)'. Then /21 = 
fii = /i3 = 1 and all other /'s equal 0, so that by (2.8), A4 = A5 = Ag = 1 
and all other A's are zeros. Hence by (2.9) and (2.10), 

^1=4, ^2=3, 13 = 2, ^4 = 3, ^5 = 4, ^6 = 6, 

h = 2, ^8 = 1, ^9 = ^10 = 3, Pi= P2 = C = ^■ 

As a result, parts (a), (b) and (d2) of Theorem 2 entail 48 words each with 
aliasing index ^; of these, 8 have length four, 32 have length five and 8 have 
length 6. Similarly, part (c) of Theorem 2 entails three words having lengths 
six, eight, and six, and each with aliasing index 1. Hence by (1.2) and (1.3), 
in this case the QC design D, which is a 2^^~^ design, has resolution 4.5 
and WLP (0,0,0,2,8,4,0,1,0,0). As seen later in Table 3, this design has 
maximum resolution and MA among all 2^^~'^ QC designs. Also, it has the 
same WLP but higher resolution than the regular 2^""^ MA design. 

Even though the trigonometric formulation keeps our derivation tractable. 
Theorem 2 is considerably more involved than its counterpart, namely, The- 
orem 1 of [16], for quarter fractions. Consequently, analytical expressions 
for the optimal choice of u and v, or equivalently, of Ai, . . . , Aio, maximizing 
the resolution or minimizing the aberration of D do not exist in easily com- 
prehensible forms. On the other hand, as Example 1 demonstrates, for any 
given Ai,...,Aio, the resolution and WLP of D can be obtained immedi- 
ately from Theorem 2. Hence, we find the best choice of the A's, with regard 
to resolution and aberration, by complete enumeration of all possible non- 
negative integer- valued Ai, . . . , Aio subject to Ai + h Aio = n, a condition 

which is evident from (2.7) and (2.8). Because of the substantial reduction 
of the problem as achieved in Theorem 2, such complete enumeration can 
be done instantaneously, for example by MATLAB, for reasonable values of 
n. The results are summarized in Table 3 and discussed in the next section 
along with their counterparts for QC designs in 2^""'"^ runs. 
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2.2. One-eighth fractions. Deletion of any one of the first four columns 
of the matrix D in Section 2.1 leads to a QC design involving 2n + 3 two-level 
factors and 2^" runs, that is, a 2(^""'"^)~^ QC design. The competing class 
of QC designs, corresponding to all possible choices of u and v, remains the 
same up to isomorphism whichever of these four columns is deleted. This 
follows by interchanging the roles of u and v and noting that 



cos 



/vr TT . \ /vr u \\ 

(M»'(3.)) 



cos I + V j = sm 

as a'u and a'v are integers. Therefore, without loss of generality, we consider 
the deletion of the first column of D. Let D^^^ denote the resulting design. 
Since the deletion of the first column of D amounts to dropping the factor 
-Fi, continuing with the notation of Section 2.1, only collections of factors 
of types X = 0000,0001,0010,0011,0100,0101,0110 or 0111 can arise now. 
For all such x, Theorem 1 again describes the numbers of words of type x 
as well as the aliasing indices and lengths of these words. Analogously to 
Theorem 2, this can be summarized as follows. 

Theorem 3. With reference to the 2(2"+^)"^ QC design D^^\ the fol- 
lowing hold: 

(a) There are 2/ p\ words each with aliasing index pi; half of these words 
have length Zi + 1 and the rest have length Z3 + 3. 

(b) There are 2/ P2 words each with aliasing index p2 and length /2 + 1. 

(c) There is one word with aliasing index 1 and length /g + 2. 

(d) In addition, 

(dl) if A5 + Ae = 0, then there are 2/(^1^2) words each with aliasing 

index £,1^,2 and length Ig + 2; 
(d2) if A5 + As > 0, then there are 2/^'^ words each with aliasing index 

^; half of these words have length Zg + 2 and the rest have length 

ho + 2. 

Theorem 3 shows that the resolution and WLP of D^^^ depend on u and 
V only through Ai, . . . , Aio and greatly simplifies the task of finding, by com- 
plete enumeration, the optimal A's maximizing the resolution or minimizing 
the aberration of D^^\ The results are summarized in Table 4 and discussed 
in the next section. 

3. Quaternary code designs in 2^"'+-'^ runs. First, consider one-sixteenth 
fraction QC designs in 2^""^^ runs as obtained by a branching technique stud- 
ied in [16] for quarter fractions. In the present context, this technique can be 
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conveniently described as follows. Let u = {uq, ui, . . . , Un)' , v = {vq, vi, . . . , Vn)' 
be (n + 1) X 1 vectors and In+i be the identity matrix of order n + 1 over 
Z4. Consider the QC given by the generator matrix [u v In+i], and let C 
be a collection of 2^""''^ codewords thereof, each of size n + 3, as given by 



(3.1) 



C = {{d'u,d'v,ao,ai, . . . , a„) : = 0, 1; ai, . . . , a„ G Z4}, 



where d = (ao,ai, . . . ,a„)', and d'u and d'v are reduced mod 4. Apply the 
Gray map (2.2) to C to get a 2^^^^ x (2n + 6) matrix D having entries ±1. 
By (3.1), the entries in the fifth and sixth columns of D correspond to the 
third entry oq in the codewords of C. Since qq = or 1, it is evident from 

(2.2) that every entry in the fifth column of D is 1, while in the sixth column 
of D half of the entries equal 1 and the remaining half —1. Delete the fifth 
column of D to get the final design matrix D^, of order 2^"+^ x (2n + 5) and 
having entries ±1. With its columns and rows identified with factors and 
runs, respectively, Dq represents a design involving 2n + 5 two-level factors 
and 2^"+^ runs. In this sense, will be called a 2(^"+^)~^ QC design. 
Evidently, the role of uq and vq in this construction is different from that of 
ui, . . . ,Un and vi,. . . ,Vn and this will be reflected in the statistical properties 
of Do- Let u = (ui, . . . , Un)' and v = {vi, . . . , Vn)'- 

We consider a trigonometric representation for the runs in Dq. Since Dq 
is obtained by deleting the fifth column of D, by (3.1) and analogously to 

(2.3) , Do has 2^" runs 



V2 



TT vr 



vr TT 



4 + 2*^^ )>cosl- + -an l,sml - + -ai; 



TT TT 



(3.2) 



TT TT 



TT TT 



,os[- + -avj,—,smy- + -a. 



TT TT \ I TT TT \ I TT TT 

cos( 4 + 2*^1 l,---,sin( - + -a„ l,cos( - + -a. 



ai, . . . ,a„ G Z4, 



which correspond to oq = 0, and another 2^" runs 
^/2 



(3.3) 



.TtT TT . ] ( TT TT . ] 

siM - + -{uo + a u) '>,cos< J + -{uo + a u) '>, 

.fvr TT , } ( TT TT , /] 1 

sin|- + -(vo + an)|,cos|- + -(uo + ai;)|,--^, 

/ TT TT \ 



f TT TT 

sm( - + -ai ),cos 



vr TT \ f TT TT 

sin( - + -«„ l,cos( - + -a„ 



fli, . . . ,a„ G Z4, 
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which correspond to ao = 1. Here a = (oi, . . . , a„)'. Denote the 2n + 5 factors 
in Do by Fi , . . . , F4 , F5 , Fii , Fi2 , . . . , Fni , Fn2 in conformity with the ordering 
in (3.2) or (3.3). 

In the spirit of Section 2, for any nonempty coUection of factors, let = 1 
if -Ffc occurs in the cohection, and otherwise, 1 < A; < 5. The coUection is 
then said to be of the type {x,x^), where x = X1X2X3X4. Thus, a typical 
collection of type {x,X5) consists of factors Fk with Xfc = 1 (1 < A; < 5), and 
also factors Fji,Fj2 {j e Si),Fj2 {j G S2) and Fji {j £ S^), where 81,82,33 
are any disjoint and possibly empty subsets of {1, ... ,n}. Such a collection 
has m + X + x^ factors, where X and m are as in Section 2, and by (1.1), 
(3.2) and (3.3), its aliasing index equals \G{x) + (—1)^^ H{x)\, where 



with (j){x;ai, . . . ,an) defined by (2.5) and (j)*{x;ai, . . . ,an) defined similarly 
replacing a'u and a'v in (2.5) by uo + a'u and vo + a'v. Any choice of 81,82, 83 
making G{x) + {—l)^^H{x) nonzero entails a word of length m + X + x^ and 
aliasing index \G{x) + {—l)^^H{x)\. Such a word is called a word of type 

{X,X5). 

Steps similar to but more elaborate than those in the Appendix may now 
be employed to develop an analog of Theorem 1 giving an account of words 
of all possible types for the 2^^""'"^)"^ QC design Dq. As hinted above, this 
has to be done separately for each possible pair uoVq. One can, thereafter, 
summarize the findings to get a counterpart of Theorem 2. However, given 
the multitude of possibilities for the pair uqVq, a tabular representation of 
these summary results is easier to comprehend than a statement in the form 
of a theorem. For any uqVq and any combination of the wordlength (wl) 
and aliasing index (ai), denote the corresponding number of words in Dq 
by N{uoVo, wl, ai)/{ai)'^. Table 1 lists all possible {wl, ai) and, for any such 
{wl, ai), shows A^(mo^O) wl, ai) for every uqVq. The derivation underlying this 
table is omitted to save space. In Table 1, h, ■ ■ ■ ,/io ai'e as in (2.9), where 
Ai, . . . , Aio continue to be given by (2.8) with reference to u = {ui, . . . , Un)' 
and V = {vi, . . . , Vn)' ■ Also, 

di = l/2<(V2)(Ai+A3+A5+A6+5i)>^ ^ ^^2({l/2)(A2+A4+A5+A6+52)>^ 



(3.4) iOo = UJIUJ2, ^1 = 1 /2<(V2)(Ai+A3+eO> ^ = i /2((l/2)(A2+A4+.2)> 
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Table 1 

Values of N{uovo,wl, at) for the 2(2"+5)-4 qq design Do 



UoVo 



wl, ai 


00 


01/03 


02 


10/30 


11/33 


12/32 


13/31 


20 


21/23 


22 


h + l,Oi 


2 


2 


2 


1 


1 


1 


1 











li+2,6i 











1 


1 


1 


1 


2 


2 


2 


J_ 1 fi^ 

t2 -r i-,fy2 


9 


1 
1 


n 
U 


9 


1 


u 


1 


9 


1 


n 
u 


h + 2, 02 





1 


2 





1 


2 


1 





1 


2 


13+3,01 


2 





2 


1 


1 


1 


1 





2 





h + 4:,di 





2 





1 


1 


1 


1 


2 





2 


k + 3, 02 


2 


1 








1 


2 


1 


2 


1 





k + 4, 02 





1 


2 


2 


1 





1 





1 


2 


«5+2,l 


1 


1 


1 














1 


1 


1 


h+3,1 











1 


1 


1 


1 











k+2,1 


1 





1 


1 





1 





1 





1 







1 








1 





1 





1 





l7+i,l 


1 





1 





1 





1 


1 





1 


l7+5,l 





1 





1 





1 








1 





Is + 2,ajo* 


4 


2 





2 





2 








2 


4 


Is + 3,LiJ0* 





2 


4 


2 





2 





4 


2 





h+2,L0* 


2 


1 





1 





1 


2 





1 


2 


^9 +3,CJ* 





1 


2 


1 


2 


1 





2 


1 





ho + 2,uj* 


2 


1 





1 


2 


1 








1 


2 







1 


2 


1 





1 


2 


2 


1 





* Entries in 


these 


rows, 


except 


the ones 


for uqvo 


= 11,13,31,33, arise 


if and only if 



As + Ae =0; 

^Entries of these rows, except the ones for uqVq = 11,13,31,33, arise if and only if 
As + Ae > 0. 

where 5i = /(uq = 1 or 3), 82 = I{vq = 1 or 3), ei = I{uoVo = 10, 12,30 or 32), 
£2 = I{uqVo = 01, 03, 21 or 23), e = ei + £2, and /(•) is the indicator function. 

Example 2. With n = 2, let n = (l,2)',u = (2, 1)' and uo'Vo = 11. Then 
/12 = /21 = 1 and all other /'s equal 0, so that by (2.8), A3 = A4 = 1 and all 
other A's are zeros. Hence, by (2.9) and (3.4), 

^1=^2 = 3, l3 = U = l, 1^ = 1^ = 2, h = ^i Is = 19 = ho = 2, 
^1=^2 = ^, Wo = 1, ^=^- 

As a result, Table 1 shows that there are 48 words each with aliasing index 
^; of these, 24 have length four and 24 have length five. In addition, there 
are three words having lengths five, five and eight, and each with aliasing 
index 1. Hence by (1.2) and (1.3), in this case the QC design Dq, which is 
a 2^-4 design, has resolution 4.5 and WLP (0,0,0,6,8,0,0,1,0,0). As seen 
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later in Table 3, this design has maximum resolution and MA among all 
2^~^ QC designs. Also, it has the same WLP but higher resolution than the 
regular 2^~^ MA design. 

We next turn to one-eighth fraction QC designs in 2^""''^ runs. Deletion of 
any one of the first four columns of Dq introduced earlier leads to a design 
involving 2n + 4 two-level factors and 2'^^'^^ runs, that is, a 

2(2n+4)-3 QC 

design. Using the same logic as in Section 2.2, without loss of generality, 
suppose the first column of Dq is deleted. Let D^^^ denote the resulting 
design. Table 2 lists all possible combinations {wl, ai) of the wordlength and 
aliasing index in Dq^^ and, for any such {wl,ai), shows N(uqVq,wI, ai) for 
every uqVq, where N(uoVo, wl,ai) is defined as above but now refers to 
In Table 2, li,l2 etc. are as in (2.9), 9i,92,ujo ai^d uj are as in (3.4), and 

A;ii = i/(Ai + A3>0), ki2 = l-kii, 
A;2i=/(Ai + A3 + A5 + A6>0), A;22 = 2-/c2i. 

Table 2 

Values of N{uovo, wl, ai) for the 2'^"'+'*'^'^ quaternary code design Dq'^ 



UoVo 



wl, ai 


00 


01/03 


02 


10 


11 


12 


13 


20 


21/23 


22 


30 


31 


32 


33 


^1 + 1,^1 


1 


1 


1 


fell 


fell 


fell 


fell 











ki2 


ki2 


ki2 


ki2 


^1+2,^1 











ki2 


ki2 


ki2 


ki2 


1 


1 


1 


fell 


fell 


fell 


fell 


12 + 1,02 


2 


1 





2 


1 





1 


2 


1 





2 


1 





1 


h + 2, 92 





1 


2 





1 


2 


1 





1 


2 





1 


2 


1 


l3 + 3,9i 


1 





1 


fell 


ki2 


fell 


ki2 





1 





ki2 


fell 


ki2 


fell 







1 





ki2 


fell 


ki2 


fell 


1 





1 


fell 


ki2 


fell 


ki2 


k + 2,1 


1 





1 


1 





1 





1 





1 


1 





1 





k + S,l 





1 








1 





1 





1 








1 





1 


Is + 2,ajo* 


2 


1 





k2i 





fc22 








1 


2 


^22 





k2i 





Is + 3,ujo* 





1 


2 


k22 





k21 





2 


1 





k21 





k22 





l9+2,U}* 


1 


1 

2 





1 

2 





1 
2 


1 





1 

2 


1 


1 
2 


1 


1 

2 





l9+3,U}* 





1 

2 


1 


1 
2 


1 


1 
2 





1 


1 

2 





1 
2 





1 
2 


1 


llO + 2,LO* 


1 


1 

2 





1 
2 


1 


1 
2 








1 

2 


1 


1 
2 





1 
2 


1 


llO+3,LO* 





1 
2 


1 


1 
2 





1 
2 


1 


1 


1 
2 





1 
2 


1 


1 
2 






*Entries in these rows, except the ones for uqvo = 11,13,31,33, arise if and only if 
A5 + A6=0; 

^Entries of these rows, except the ones for uqVo = 11,13,31,33, arise if and only if 
As + Afi > 0. 
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Table 3 

One-sixteenth fraction QC designs with maximum resolution and MA 



Regular 



Design 




QC design with mciximum resolution and MA 




MA design 


28-4 


A 


= 0011000000, 7? = 4, A = (14, 0, 0, 0, 1) 


R 


= 4, A same 


29-4 


A 


= 0011000000, novo = 11,7? = 4.5, A = (6, 8, 0, 0, 1, 0) 


R 


= 4, A same 


2IO-4 


A 


= 0001110000, R = 4.5, yl = (2,8,4,0,1,0,0) 


R 


= 4, A same 


2II-4 


A 


= 0001110000, novo = 12, R = 5.5, A = (0, 6, 6, 2, 1, 0, 0, 0) 


R 


= 5, A same 


2I2-4 


A 


= 0011110000, J? = 6.5, yl = (0,0,12,0,3,0,0,0,0) 


R 


= 6, A same 


213-4 


A 


= 0011110000, UQVQ = 22, 7? = 6.5, A = (0, 0, 4, 8, 3, 0, 0, 0, 0, 0) 


R 


= 6, A same 


214-4 


A 


= 1011110000, J? = 6.5, 4= (0,0,2,8,3,0,2,0,0,0,0) 


R 


= 7, A better 



As illustrated in Example 2, Tables 1 and 2 readily yield, in their re- 
spective contexts, the resolution and WLP of a QC design for any given 
Ai,...,Aio and uqVq. Hence, we find the best choice of the A's and uqVq, 
with regard to resolution and aberration, by complete enumeration of all 
possibilities. Again, because of the significant reduction achieved in Tables 
1 and 2 by theoretical means, such complete enumeration can be done in- 
stantaneously, for example by MATLAB, for reasonable values of n. The 
results are summarized in Tables 3 and 4 for one-sixteenth and one-eighth 
fractions, respectively. 

A brief discussion on Tables 3 and 4, showing one-sixteenth and one- 
eighth fraction QC designs with maximum resolution and MA is in order. 
The 14 QC designs shown in these tables are optimal, among all comparable 
QC designs, under both criteria. All these 14 designs have resolution four or 
higher, and for each, the resolution R and A = (^44, ^5, . . .) are shown. For 
ease in comparison, we also show R and comment on A for the corresponding 
regular MA designs, as obtained by Chen and Wu [6]. Out of the 14 optimal 
QC designs in Tables 3 and 4, there are two, that is, the first design in 

Table 4 

One-eighth fraction QC designs with maximum resolution and MA 



Regular 



Design 




QC design with maximum resolution and MA 


MA design 


27-3 


A = 


0011000000, i? = 4, A = (7, 0,0,0) 


R 


= i, A same 


28-3 


A = 


0011000000, tiouo = ll, i? = 4.5, A = (3,4,0,0,0) 


R 


= 4, A same 


29-3 


A = 


0010110000, 7? = 4.5, A = (1,4,2,0,0,0) 


R 


= 4, A same 


2IO-3 


A = 


0010110000, liouo =21, i? = 5.5, A = (0, 3, 3, 1,0,0,0) 


R 


= 5, A same 


2II-3 


A = 


0011110000, R = 6.5, A = (0,0,6,0,1,0,0,0) 


R 


= 6, A same 


2I2-3 


A = 


0011110000, ^01)0 = 12, i? = 6.75, A = (0,0,2,4,1,0,0,0,0) 


R 


= 6, A same 


2I3-3 


A = 


0021110000, R = 7.75, A = (0,0,0,4,3,0,0,0,0,0) 


R 


= 7, A same 
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either table, which have the same R and A as the corresponding regular 
MA designs. It can be seen that these two designs involve only full words 
and hence are themselves regular. A comparison with Table 6 of Sun, Li 
and Ye [18] shows that these two designs have maximum resolution and MA 
in their sense among all designs of the same size. Eleven of the remaining 
twelve optimal QC designs in our Tables 3 and 4 have the same WLP but 
higher resolution than the corresponding regular MA designs. Only the 2^^~^ 
optimal QC design turns out to be worse than the regular MA design. In 
both Tables 3 and 4, A stands for the 10-tuple A1A2 • • • Aiq. 

Indeed, the theoretical results reported in Theorems 2, 3 and Tables 1, 2 
readily allow extension of Tables 3 and 4 beyond the ranges considered 
here, if the situation so demands. For instance, from Table 2, one can check 
that the 2^^^'^ QC design with maximum resolution and MA is given by A = 
0020220000, uovo = 20. This design has R = 8.875 and A = (0, 0, 0, 0, 1, 4, 2, 0, 
0,0,0,0,0), while the corresponding regular MA design has the same A but 
R = 8. 

4. Results on projectivity. The following results give upper bounds on 
the projectivity of the one-sixteenth fraction QC designs D and Dq intro- 
duced in Sections 2 and 3, respectively. 

Theorem 4. The projectivity p of the 2(2"+"^)"^ QC design D satisfies 
(i) P < |n + 1, if n = 0mod3, (ii) p < |(n — j) -|- 3, if n = j mod 3, with j = 1 
or 2. 

Proof. We prove only (ii). The proof of (i) is similar. Let n = 3t + j, 
where t is an integer and j = 1 or 2. We need to show that p < 4t + 3. If 
p> 4t + 4, then all full words in D have length at least 4t + 5, so that by 
Theorem 2(c) and (2.9), 

15 + 2 = 2(Ai + A3 + As + Ae) + 2 > 4t + 5, 

k + 2 = 2(A2 + A4 + As + Ae) + 2 > 4t + 5, 

/7 + 4 = 2(Ai + A2 + A3 + A4) + 4 > 4t + 5, 

that is, invoking the integrality of Ai, . . . , As, 

Ai + A3 + As -h Ae > 2t + 2, A2 + A4 + As -h Ae > 2t + 2, 

Ai + A2 + A3-hA4>2t + l. 

Adding the last three inequalities, 2(Ai + • • • + Ag) > 6t + 5, that is, n > 

Ai H \- XQ>3t + 3, again using the integrality of Ai, . . . , Ag, and we reach 

a contradiction. □ 



16 



ZHANG, PHOA, MUKERJEE AND XU 



From Table 1, now observe that the 2^^""*"^^ ^ QC design Dq has at 
least three full words. These have lengths (a) + 2, /g + 2, /y + 4 if uqVq = 
00, 02, 20, 22, (b) k + 2, + 3, /t + 5 if uqVq = 01, 03, 21, 23, (c) k + 3, + 2, 
/7 + 5 if uqVo = 10, 12, 30, 32 and (d) ^5 + 3, Zg + 3, ^7 + 4 if uqVo = H, 13, 31, 33. 
Hence, arguments similar to but more elaborate than those in Theorem 4 
lead to the following result. 

Theorem 5. The projectivity p of the 2(2'^+5)-4 qq design Dq satisfies 
(i) p < |n + 2 if n = 0mod3, (ii) p < |(n — 1) + 3, if n = lmod3, (iii) 
p< |(n-2) +4, i/n = 2mod3. 

Table 5 shows the projectivities of the one-sixteenth fraction QC designs 
reported in Table 3. It is easily seen that these designs attain the upper 
bounds on projectivity as shown in Theorems 4 or 5. Hence, in addition 
to having maximum resolution and MA, they have maximum projectivity 
among all comparable QC designs. Indeed, for 8 < g' < 12, the 2'^"^ QC 
designs in Table 3 have projectivity q — b which is the highest among all 
designs of the same size. This holds because, otherwise, one would get a 2*?"^ 
design with projectivity <? — 4, that is, an orthogonal array OA (2''"^, q,2,q — 
4) of index unity, which is nonexistent; see [3]. Table 5 also shows that the 
use of QC designs leads to gain in projectivity over regular MA designs for 
9 < g < 14. 

We next consider one-eighth fractions and show in Table 6 the projectiv- 
ities of the QC designs reported in Table 4. For 7 < g < 11, the 2''"^ QC 
designs in Table 4 are seen to have projectivity q — 4 which is the highest 
among all designs of the same size. This follows as in the last paragraph 
using a nonexistence result in [3] on orthogonal arrays of index unity. Also, 
for q = 12 and 13, the QC designs in Table 4 were computationally verified 
to have maximum projectivity at least among all comparable QC designs. 
Incidentally, for one-eighth fraction QC designs, it is hard to develop analogs 
of Theorems 4 and 5 as there is only one guaranteed full word (cf. Theo- 
rem 3 and Table 2) but the computational study of projectivity remains 



Table 5 

Projectivities of the one-sixteenth fraction QC designs in Table 3 



Design 2^-* 2""^ 


210-4 


211-4 


212-4 


213-4 


214-4 


Projectivity of the 3 4 


5 


6 


7 


7 


7 


QC design in Table 3 












Projectivity of the 3 3 


3 


4 


5 


5 


6 



regular MA design 
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Table 6 



Projectivities 


of the 


one-eighth 


fraction 


QC desi 


gns in Table 4 




Design 


2T-3 


28-3 


29-3 


2IO-3 


2I1 — 3 2^-^"^ 


213-3 


Projectivity of the 


3 


4 


5 


6 


7 7 


7 


QC design in Table 4 














Projectivity of the 


3 


3 


3 


4 


5 5 


6 



regular MA design 



manageable with a moderate number of factors. Table 6 also shows the pro- 
jectivity of regular MA designs and the gains via the use of QC designs, for 
8 < q < 13, are evident. 

The foregoing discussion reveals that, unlike what often happens with 
quarter fraction QC designs [16], maximum projectivity is not in conflict 
with maximum resolution or MA in our setup at least over the range covered 
by Tables 3-6, where the same QC design turns out to be optimal, among 
all comparable QC designs, with regard to all the three criteria. 

5. Summary and future work. In the present paper, a trigonometric ap- 
proach was developed to obtain theoretical results on QC designs with focus 
on one-eighth and one-sixteenth fractions of two-level factorials. It was seen 
that optimal QC designs often have larger resolution and projectivity than 
comparable regular designs. In addition, some of these designs were found 
to have maximum projectivity among all designs. 

Before concluding, we indicate a few open issues. It should be possible 
to use the trigonometric approach to obtain further theoretical results on 
the projectivity of QC designs. This calls for examining the existence of 
solutions to certain trigonometric equations. For instance, by (2.3), the first 
four columns of the 2(^'"+^)~^ QC design D in Section 2 contain a full 2^ 
factorial if and only if for every yi, ... ,1/4 in { — 1,1}, the equations 




admit a solution for a = (ai,...,an)' in Z4. A study of equations of this 
kind, however, branches out into too many cases, depending on u and v, 
compared to the derivation of results on wordlength and aliasing index as 
done here. 

It would also be of interest to investigate how the trigonometric approach 
can be implemented for even more highly fractionated QC designs. The 
trigonometric formulation as well as the mathematical tools are expected 
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to be essentially same as the ones here and the main difficulty will lie in 
handling the multitude of cases that such an effort will involve. A related 
issue concerns the development of a complementary design theory for QC 
designs in the spirit of [5, 21] and with respect to an appropriately defined 
reference set. 

Some kind of symbolic computation may help in addressing the open 
problems mentioned above. We hope that the present endeavor will generate 
further interest in these and related issues. 

APPENDIX: PROOF OF THEOREM 1 FOR X = 0101 

The proof will be worked out through a sequence of lemmas. We concen- 
trate on V{x) and begin by giving an expression for V'(ai) • • • i^n) in (2-6). 
Recall that in (2.6), 51,52,53 are disjoint subsets of {l,...,n} and that 
m = 2n\ + 77,2 + with = #5/;. For 1 < A; < 3, let denote the sum over 
2*^* terms corresponding to the 2*^* subsets Wk of 5/;, and for any such sub- 
set Wk, write Wk = Sk\Wk,Wk = #Wk,Wk = #Wk. Thus if 5i = {2, 3}, then 
Si denotes the sum over 2^ terms corresponding to Wi = empty set, {2}, 
{3} and {2, 3}. For any given subsets VFi, VF2, VF3 of 5i, 52, 53, we also write 
S(1),S(1\S(23^ and S^^a) denote sums over j G Wi,j G Wi,j G W2 U Ty3 
and j G W2 U W^, respectively. Similarly, S^^^ denotes sum over j £ 84^, where 
54 = {1, . . . , ?i} \ (5i U 52 U 53). Let i = \/— 1. Then the following lemma is 
not hard to obtain using elementary facts such as sin y cos y = ^ sin 2y, cos y = 
\ {e^y + e-'y), sin y = ^ {e^y - e'^y) . 

Lemma A.l. 

1p{ai,...,an) = 2mjli+n3 '^1^2^3M{wi , -W2 , W3) 

X exp|^(2S(i)a, - 2t('^a, + E^^^^a, - S(23)a,.)|, 
where M{wi,W2,W3) = (-l)'^i+'^3 exp{^(m - 4u)i - 2tt;2 - 2tt;3)}. 

Let g = {gi,...,gn)' ,h = {hi,...,hn)\ where 

(A.l) gj = Uj + Vj{m.odA)^ hj = uj — Vj{mod4:), 1 < j < ra. 

Using the same elementary facts that led to Lemma A.l, we now note 
that 

/ vr TT , \ / vr TT , 

cos — I — a u cos — I — a v 

V4 2 y V4 2 



(A.2) 



exp|y(l + a'g)| -hexpjya /i| 

I ZTT 1 1 ITT 1 
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Throughout the rest of the Appendix, including the lemmas below, we con- 
sider X = 0101. Then X = xi + X2 + + = 2, and hence (2.5), (A. 2) and 
Lemma A.l yield 

(A.3) 

4 



^(/)^(x;ai,...,an), 



X 
k=l 



where, with gQ = 1, 

(l)Y{x;ai, ...,an) 



(A.4) = exp 



ITT 

y 



{go + S(i)(<7, + 2)aj + J^^'^gj " 2)a 



'J 



and, for k = 2,3, 4, (/'^(x; ai, . . . , 0^) are analogous to i^Y {^'i di, - ■ ■ , a^), with 
(50,51, • • ■,an) in the latter replaced by (0,/ii, . . . ,/in), (0,-/ii,. . .,-hn) and 
(—go, —51, • • • , —9n), respectively. The superscript W here indicates the de- 
pendence of each (/)^(x; ai, . . . , a„) on Wi,W2,W3 via the sums S^^^S^^^ 
5^(23)^ 5^(23) _ Writing 

^.(-) = 2(l/2)(n.+2U^n.+n3 ^1^^^3M(u-T,^2,^3)^.^(x), 

where 

3 3 

(A.6) Vj^{x)=^---J2^k'{x;ai,...,an), 1 < k < A, 

ai=0 a„=0 

the following lemma is immediate from (2.4) and (A.3). 

Lemma A. 2. V{x) = Vi{x) + V2ix) + Vsix) + V^ix). 

Some more notation will help in presenting the subsequent lemmas. With 
reference to the sets A^^ in (2.7) and the gj and hj in (A.l), let 

A(^) = Aio U A12 U A30 U A32, A(2) = Aoi U Ao3 U A21 U A23, 

(A.7) 

a = aWua(2), 

= {j : 1 < J < n,gj = k}, A^ = {j : 1 < j < n,hj = k}, 

(A.8) 



(A.ll) 
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Then by (2.7), (2.8) and (A.l), 

/3i = #A(i) = /lo + /i2 + /so + /32 = Ai + A3, 
(A.9) h = #A(2) = /oi + /o3 + /21 + /23 = A2 + A4, 
/3 = #A = /3i + /32 = Ai + A2 + A3 + A4, 
Ag = Aoo U Ai3 U A22 U A31, A[J = Aoo U An U A22 U A33, 
(A.IO) Af = A02UA11UA20UA33, A^ = Ao2UAi3UA2oUA3i, 
Af U Af = A^ U A^' = A, 

/| = A5 + A7 + A8, 

fh = X^ + h + \^, ff + fl = f^ + f^ = l3. 
We also write 

(A. 12) (31= fu + /so - /lO - /32 , f3*2 = /03 + /2I - /Ol - /23 , 

(A.13) Mi = #(53nA«), ^2 = #(S3nA(2)), 

(A.14) n{g,k) = #{S3nAl), n(/i, A:) = #(^3 n A^!), < < 3. 

Lemma A. 3. (i) Vi{x) = V4{x) = 0, unless 
(A.15) Si = Af, 82^83 = AfuAl, 54 = Ag. 

(ii) V2(x) = V3(x) = 0, unless 
(A.16) 5i = A^, 52U53 = A?UA^, 54 = A|}. 

Proof. We prove (i). The proof of (ii) is similar. Since for any integer k, 

3 



y~^exp[ —ks ) = 4 ifA; = 0mod4, 



■5=0 

(A.17) 

= otherwise, 

it follows from (A. 4) and (A. 6) that (x) vanishes, for every Wi (cS'i), 
W2 {C82) and W3 {C83), and hence by (A.5) ^1(2;) =0, if either (a) 81 is 
nonempty and gj 7^ 2 for some j £ 81, or (b) ^2 U ^3 is nonempty and gj ^1,3 
for some j £ 82U 83, or (c) ^4 is nonempty and gj 7^ for some j £ 84. Thus 
Vi{x) = 0, unless 81 C Ag, 5*2 U C Af U A3 and 5*4 C Aq. These conditions 
are equivalent to those in (A.15) because 51,52,53,54 form a partition of 
{1, . . . , n} as Aq, Aj', A2, A3 do. The same arguments apply to V4(x). Hence, 
(i) follows. □ 
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Lemma A. 4. (i) // (A. 15) holds then 
Vi{x) + Vi{x) 
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2(l/2)/3+ljn3 



ITT 



exp^-(2 + /|-/f: 



ITT 



+ (-irexp^--(2 + /|-/f) 



(ii) If (A. 16) holds, then 

V2{x) + V3{x) 



_l\n{h,l) 



2(l/2)/3+l^n3 



+ (-irexp^-^(/^-A'^: 



Proof. We prove (i), the proof of (ii) being similar. Let (A. 15) hold. 
Then by (A. 4), noting that qq = 1, 



(pY{x;ai, . . . ,a„) 



(A.18) 



exp 



^{1 + S(^3)(,, + I)a, + S(2=')(ff,-I)a,} 



ITT , 



■ exp<j + 2i;*aj + 2J:**aj] 



S* and S** being sums over j G (Wa U W3) n A? and j G (W2 U VF3) n Af , 
respectively. By (A. 6), (A. 17) and (A.18), V^{x) = unless the ranges of 
S* and S** are both empty. From this, invoking the second equation in 
(A. 15), a little reflection shows that Vi^{x) = unless 



(A.19) 



W2 = S2nAl 



W3 



W2 = S2nAl W'3 = 53nAf. 



Given ^2,53, (A.19) determines W2,W3,W2 and W3 uniquely. Moreover, if 
(A.19) holds then W3 = n{g, 1) by (A.14), (x) = 4"exp(f ) by (A.6) and 
(A.18), and 



m — iwi — 2w2 — 

= 2ni — Awi + W2 + w-i — W2 — W3 

= 2ni - Awi + #{(S2 u 53) n Af } - #{(^2 U S3) n A?} 

= 2ni - 4tDi + #Af - #A? = 2ni - Awi + /| - /f 
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by (A. 8), (A. 15) and the facts that m = 2ni + n2 + n^-,nk = Wk + Wk {k = 
2,3). If we summarize the above and recah the definition of M{w\,W2,W'i) 
from Lemma A.l, then from (A. 5), we get 

^ 4"exp(z7r/2) , .u„+n(g,l) 



(A.20) 

ITT 



xexp<|^(2ni-4«;i + /|-/f)|. 



Since m + 2 = 2ni + /f + /| + 2 = 2ni + /? + 2 by (A.ll) and (A.15), and 

Ei(-l)-^exp|^(2ni-4«)i)|=(2i)"i 

as one can verify after a httle algebra, (A.20) yields 

V, (X) = j^jj^^ expj -(2 + /| - /f ) | . 

Similarly, it can be shown that under (A.15), 

(A-22) ^4(x) = ^y^exp|--(2 + /I - /f )|. 

By (A. 14) and (A.15), n3 = n{g,l) +n{g,3). Hence, (i) follows from (A.21) 
and (A.22). □ 

Lemma A. 5. (i) (-l)"(^'i)-"(ff'i) = (-l)/^2 . (ii) If either (A.15) or (A. 16) 
holds then ^3 = /^i + /i2 • 

Proof, (i) Let /xq = #{83 n (Aqi U A23)}. From (2.7), (A.l) and (A.8), 

Af = Aoi U Aio U A23 U A32, A^ = Ao3 U Aio U A21 U A32. 

Hence by (A.7), (A. 13) and (A.14), n{h,l) - n{g,l) = ^2 - and (i) 
follows. 

(ii) If either (A.15) or (A. 16) holds, then by (A.7) and (A. 10), ^3 C A 
[=A(i) U A(2)]. Now (ii) is immediate from (A.13). □ 

Lemma A. 6. (i) IfX5 + \Q = 0, then (A.15) and (A. 16) become identical 
and both reduce to 

(A.23) 5i = Ao2UA2o, 82^83 = A, 54 = AooUA22. 

(ii) If X5 + Xq> 0, then (A.15) and (A. 16) cannot hold simultaneously. 
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Proof. If A5 + Ae = 0, then by (2.7) and (2.8), the sets An, A13, A31 
and A33 are empty, so that by (A. 10) , A2 = A2 = A02 U A20 and Aq = 
Aq = Aqo U A22- Hence, (i) fohows recahing the last identity in (A. 10). On 
the other hand, if A5 + Ag > 0, then at least one of An, A13, A31 and A33 is 
nonempty. Therefore, Ag 7^ Aj and Aq ^ Aq and (ii) follows. □ 

Lemma A. 7. Let X5 + Xq = 0. Then: 

(i) V{x) = 0, if (A. 23) does not hold. 

(ii) vlx) = (-l)"(3-i)Q^(a;)g2(a;), if (A. 23) holds, where for k = 1,2, 



Qk{x) 



exp^^(l + C 



(A.24) 

+ (-irexp|-^(l+/3^; 

Proof. Part (i) is evident from Lemmas A. 2, A. 3 and A.6(i). To prove 
(ii), let (A.23) hold and write ai = f (2 + /| - /f),a2 = i{f^ - f^). Now, 
by Lemmas A. 2, A. 4, A. 5 and A.6(i), 

/_i\n(3,l) ... 

^(^) = 2(i/2)/3+i^n3 i^'"^ + (-l)'^^+'^^e— ^ + (-l)^^e-^ + (-l)'^^e— n 
(A.25) = -^^^^^^^{e*K+"^)/2 + (_i)Mi,-K+«2)/2} 

But i(ai + 02) = f (2 + /I - /f + f^ - ) = f (1 + PI), as one can verify 
from (2.7), (A.l), (A. 8) and (A. 12), on simplification. Similarly, ^(ai — 02) = 
^(1 + /32)- Hence, (ii) follows from (A.25) using Lemma A.5(ii) and the fact 
that /3 = /3i + /32; cf. (A.9). □ 

Lemma A. 8. (i) both f3i and (32 o.re odd, then the number of choices, 
say a, of Si, S2, which meet (A.23) and keep both /ii and fi2 odd equals 
2^-2. 

(ii) The number of choices of Si, S2, S-^ meeting (A. 16) is 2^. 

(iii) // /3 (>0) is even, then ao = ai = 2^~^ , where ctq and ai denote the 
numbers of choices of Si, S2, S3 which meet (A. 16) and keep even and 
odd, respectively. 

Proof. We prove only (i). Proofs of (ii), (iii) are similar. For a given QC 
design, (A.23) determines Si uniquely and fixes S'2US'3 at A [=A(^) U A^^), 
by (A. 7)]. Thus, by (A. 13), a equals the number of ways in which one can 
choose an odd number of elements from each of A^^^ and A^^^. Hence, a = 



24 



ZHANG, PHOA, MUKERJEE AND XU 



(2/3i-i)(2fe-i) = 2^- 
hence positive. □ 



as /3i(=#A(i)) and ^2(=#A(2)) are both odd and 



Proof of Theorem 1 for x = 0101. The proof is given separately for 
three cases corresponding to (il), (i2) and (i3) of Theorem 1. Recah that 
any choice of Si, S2, S3 which makes V{x) nonzero entails a word of type x, 
length m + X and aliasing index 

Case 1. Let A5 + Ae = 0. Then Lemma A. 7 is applicable. Thus if V{x) ^ 
then (A.23) holds, so that m + X = 2(/o2 + /20) + #A + 2 = + 2, by (2.7)- 
(2.9) and (A. 9). Hence, each word of type x has length + 2. It remains to 
show that there are 1/(^1^2) words of type x, each with aliasing index ^1^2- 
To that effect, from (A. 24) note that for any S*!, 6*2,53 meeting (A.23), 



(A.26) 

where for k ■ 



\V{x)\ = \Qiix)\\Q2{x)\, 



1,2, 



|Q,(rr)|=2-(V2)(A-i) 


cos< 




^ 2-(l/2)(/3fc-l) 


sin< 





if /ifc is even 



if fik is odd. 



First, suppose both f5i and P2 are odd. Then by (A. 9) and (A. 12), PI and 
are also both odd. For A; = 1, 2, if = 1 or 5 (modS), then |(5fc(2;)| equals 
or 2~(^/^)('^'=~^) depending on whether /i^ is even or odd, while if = 3 or 7 
(mod 8) , then the roles of even and odd fj-k are switched. Hence if both fS^ and 
/3| equal 1 or 5 (modS), then by (A.26), a choice of 5i, 5*2, S3 meeting (A.23) 
yields a nonzero V{x) and hence leads to a word if and only if /ii and fi2 are 
both odd. By (2.10), (A. 9) and Lemma A.8(i), the number of such choices is 
2/5-2 [=1/(^2^2^] and, for any such choice, \V{x)\ = 2-(V2)(/3i+/32-2) [=^^^2]. 
It is easily seen that the same holds if both (3* and /Sg equal 3 or 7 (mod 8) 
or if one of them equals 1 or 5 (mod 8) and the other 3 or 7 (mod 8). This 
settles Case 1 when /3i and (32 are both odd. Similar arguments work when 
they are both even or one of them is odd and the other even. 

Case 2. Let A5 + Ae > and Ai + A2 + A3 + A4 = 0. By (A.9)-(A.ll), 



here /^ = 0,/f = /| = /f 



= 0, and Af U A| and Af- U A^ are empty 



sets. Hence, by (A. 14), 71,3 = n{g, 1) = n{h, 1) = 0, for any 5i, 5*2, 63 meeting 
(A. 15) or (A. 16). Therefore, Lemma A. 4 yields Vi{x) + V4{x) = under 
(A. 15), and V2{x) + V3(x) = 1 under (A. 16). Since Lemma A.6(ii) rules out 
simultaneous occurrence of (A. 15) and (A. 16), Lemmas A. 2 and A. 3 show 
that V{x) equals 1 if (A. 16) holds, and otherwise. Moreover, as A^ U A3 
is empty, (A. 16) determines 5i,52,S'3 uniquely and, under (A. 16), m + X = 
2/2 + 2 = lio + 2, using (2.9), (A. 11). Thus, in this case, there is one word 
of type X, with aliasing index 1 and length Iiq + 2. 
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Case 3. Let A5 + Ae > and Ai + A2 + A3 + A4 > 0. Again, Lemma A.6(ii) 
precludes coincidence of (A. 15) and (A. 16). Thus, by Lemmas A. 2 and A. 3, 
words of type x can arise in two mutually exclusive and exhaustive ways: 
(a) from 5i, 52,^3 meeting (A. 15), (b) from Si, 5*2, 5*3 meeting (A. 16). First, 
consider (b). For any 5i, 52,53 meeting (A. 16), by Lemmas A. 2, A. 3 and 
A.4(ii), 



\V{x)\ 



-(l/2)/3 



-(l/2)/3 



COs|j(/3"-/f) 

sin|^(/3^-/f) 



if 712, is even 



if n3 is odd. 



Thus for odd ft^ — /(', irrespective of whether is even or odd, |V^(2;)| = 
2-(i/2)(/3+i)^ On the other hand, if fl^ - = or 4 (mod 8), then \V{x)\ 
equals 2"^^/^)^ or according as whether 77,3 is even or odd, while if — = 
2 or 6 (mods), then the roles of even and odd ns are switched. Since /3 
[=fi + /3, by (A. 11)] and fi^ — fi are either both even or both odd, from 
(2.10), (A. 9) and Lemma A.8(ii), (iii), it follows that, irrespective of whether 
/s ~ fi is odd or even, (b) yields 1/(2^^) words of type x, each with aliasing 
index S^. By (2.9), (A. 9), (A. 11) and (A. 16), each such word has length 
m + X = 2f^ + + f^ + 2 = ho + 2. Similarly, (a) yields another 1/(2^^) 
words of type x, each with aliasing index and length /g + 2. Hence, the 
conclusion of Theorem 1 for x = 0101 follows in this case. □ 



Acknowledgments. We thank the referees for their very helpful suggestions. 



REFERENCES 

[1] Box, G. and Tyssedal, J. (1996). Projective properties of certain orthogonal arrays. 

Biometnka 83 950-955. MR1440059 
[2] Box, G. E. P. and Hunter, J. S. (1961). The 2*"'' fractional factorial designs. 

Technometrics 3 311-351. 449-458. 
[3] Bush, K. A. (1952). Orthogonal arrays of index unity. Ann. Math. Statist. 23 426- 

434. MR0049146 

[4] Butler, N. A. (2003). Minimum aberration construction results for nonregular two- 
level fractional factorial designs. Biometnka 90 891-898. MR2024764 

[5] Chen, H. and Hedayat, A. S. (1996). 2""' designs with weak minimum aberration. 
Ann. Statist. 24 2536-2548. MR1425966 

[6] Chen, J. H. and Wu, C. F. J. (1991). Some results on s""*" fractional factorial 
designs with minimum aberration or optimal moments. Ann. Statist. 19 1028- 
1041. MR1105859 

[7] Cheng, C.-S., Deng, L.-Y. and Tang, B. (2002). Generalized minimum aberration 
and design efficiency for nonregular fractional factorial designs. Statist. Sinica 
12 991-1000. MR1947057 

[8] Cheng, S.-W., Li, W. and Ye, K. Q. (2004). Blocked nonregular two-level factorial 
designs. Technometrics 46 269-279. MR2082497 



26 



ZHANG, PHOA, MUKERJEE AND XU 



[9] Deng, L.-Y. and Tang, B. (1999). Generalized resolution and minimum aberra- 
tion criteria for Plackett-Burman and other nonregular factorial designs. Statist. 
Smica 9 1071-1082. MR1744824 

[10] Fries, A. and Hunter, W. G. (1980). Minimum aberration 2*'"*' designs. Techno- 
metrics 22 601-608. MR0596803 

[11] Hammons, A. R. Jr., Kumar, P. V., Calderbank, A. R., Sloane, N. J. A. and 
Sole, P. (1994). The Z4-linearity of Kerdock, Preparata, Goethals, and related 
codes. IEEE Trans. Inform. Theory 40 301-319. MR1294046 

[12] Li, Y., Deng, L.-Y. and Tang, B. (2004). Design catalog based on minimum G- 
aberration. J. Statist. Plann. Inference 124 219-230. MR2066236 

[13] Ma, C.-X. and Fang, K.-T. (2001). A note on generalized aberration in factorial 
designs. Metrika 53 85-93 (electronic). MR1836867 

[14] MuKERjEE, R. and Wu, C. F. J. (2006). A Modern Theory of Factorial Designs. 
Springer, New York. MR2230487 

[15] Phoa, F. K. H. (2009). Analysis and Construction of Nonregular Fractional Factorial 
Designs. Ph.D. thesis, Univ. California, Los Angeles. MR2712984 

[16] Phoa, F. K. H. and Xu, H. (2009). Quarter-fraction factorial designs constructed 
via quaternary codes. Ann. Statist. 37 2561-2581. MR2543703 

[17] Stufken, ,L and Tang, B. (2007). Complete enumeration of two-level orthogonal 
arrays of strength d with rf 4- 2 constraints. Ann. Statist. 35 793-814. MR2336869 

[18] Sun, D. X., Ll, W. and Ye, K. Q. (2011). An algorithm for sequentially con- 
structing non-isomorphic orthogonal designs and its applications. Unpublished 
manuscript. 

[19] Tang, B. (2001). Theory of J-characteristics for fractional factorial designs and 
projection justification of minimum G'2-aberration. Biometrika 88 401-407. 
MR1844840 

[20] Tang, B. and Deng, L.-Y. (1999). Minimum G2-aberration for nonregular fractional 
factorial designs. Ann. Statist. 27 1914-1926. MR1765622 

[21] Tang, B. and Wu, C. F. J. (1996). Characterization of minimum aberration 2"~'° 
designs in terms of their complementary designs. Ann. Statist. 24 2549-2559. 
MR1425967 

[22] Wu, C. F. J. and Hamada, M. (2000). Experiments: Planning, Analysis, and Pa- 
rameter Design Optimization. Wiley, New York. MR1780411 

[23] Xu, H. (2005). Some nonregular designs from the Nordstrom-Robinson code and 
their statistical properties. Biometrika 92 385-397. MR2201366 

[24] Xu, H. and Deng, L.-Y. (2005). Moment aberration projection for nonregular frac- 
tional factorial designs. Technometrics 47 121-131. MR2188074 

[25] Xu, H., Phoa, F. K. H. and Wong, W. K. (2009). Recent developments in non- 
regular fractional factorial designs. Statist. Surv. 3 18-46. MR2520978 

[26] Xu, H. and WONG, A. (2007). Two- level nonregular designs from quaternary linear 
codes. Statist. Sinica 17 1191-1213. MR2397390 

[27] Xu, H. and Wu, C. F. J. (2001). Generalized minimum aberration for asymmetrical 
fractional factorial designs. Ann. Statist. 29 1066-1077. MR1869240 

[28] Ye, K. Q. (2003). Indicator function and its application in two-level factorial designs. 
Ann. Statist. 31 984-994. MR1994738 



QUATERNARY CODE DESIGNS 



27 



R. Zhang 

Klas and School of Mathematics 



F. K. H. Phoa 

Institute of Statistical Science 

ACADEMIA SiNICA 

Taipei 11529 
Taiwan 

E-MAIL: frodphoalSistat.sinica.edu.tw 



AND Statistics 
Northeast Normal University 
Changchun 130 024 
China 



and 

LPMC AND School 

OF Mathematical Sciences 
Nankai University 
Tianjin 300 071 
China 

E-MAIL: rczhang@ncnu.edu.cn 
zhrch@nankai.edu.cn 

R. Mukerjee H. Xu 

Indian Institute of Management Calcutta Department of Statistics 

JOKA, Diamond Harbour Road University of California 

Kolkata 700 104 Los Angeles, California 90095-1-554 

India USA 

E-MAIL: rmukl@hotmail.com E-MAIL: hqxu@stat.ucla.cdu 



